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Recent angle-resolved photoemission data for overdoped Bi2212 are explained. Of the peak-dip- 
hump structure, the peak corresponds the ^=0 component of a hole condensate which appears at 
Tc. The fluctuating part of this same condensate produces the hump. The bilayer splitting is large 
enough to produce a bonding hole and an electron antibonding quasiparticle Fermi surface. Smaller 
bilayer splittings observed in some experiments reflect the interaction of the peak structure with 
quasiparticle states near, but not at, the Fermi surface. 
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The existence of a bilayer splitting remains a contro- 
versial issue. There are now two separate claims^ to have 
observed via angle-resolved photoemission spectroscopy 
(ARPES) such a splitting in modest to heavily overdoped 
Bi2212 while others^ have concluded that only a very 
small such splitting is possible for optimal or underdoped 
samples of the same material. The question of a possible 
bilayer splitting cannot be separated from the recent de- 
bate about the nature of the Fermi surface. Until recently 
the widely accepted Fermi surface was hole like^, how- 
ever there have been several reports^ that, for different 
ARPES photon energies, an electronlike Fermi surface is 
seen. Very recently for a slightly overdoped Bi2212 it has 
been claimed^ that both electron and hole Fermi surfaces 
exist and that this reflects a modestly large bilayer split- 
ting. These most recent experiments would seem to be 
in conflict with the earlier observation of a bilayer split- 
ting. In the experiments^ of Feng et ai, on a heavily 
overdoped sample the antibonding band was observed to 
cross the Fermi surface near [tt, 0] while the existence of 
an electronlike antibonding Fermi surface for a smaller 
hole doping is at odds with this observation. 

In the picture developed here, as proposed^ by Bog- 
danov et al.^ there exists a bilayer splitting which is suffi- 
ciently large that the bonding band is hole, while the an- 
tibonding band is electronlike. That part of the hole like 
antibonding band seen^ by Feng et a/., which closes near 
[tt, 0], reflects, in fact, a near Fermi surface "shadow" res- 
onance associated with an antibonding band which itself 
never crosses the Fermi surface. 

For a bilayer the usual t — J model becomes. 
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where c]^^ = (l-n. 

Acra electron of spin cr site i and plane a. The 

in-plane interactions are t and J and near- neighbor. The 
interplane exchange is vertical, while the t^ leads to 
a single-particle bilayer splitting {t^/ A){co^kx — cos ky)'^. 

The present formalism signals a departure from the 
usual slave-boson method^'^. Used here is a formalism 
based upon the Jordan- Wigner (JW) transformation^. 
The system is mapped to a one dimensional path. A 



single spin fermion, fj creates an up spin at path site 
i when it acts on the down spin ferromagnetic vacuum 
\)-N/2' Here the total spin Sz — —N/2 and N is the 
(even) number of sites. A hole is created hy h\. A hard- 
core constraint = f} h\h^ < 1 is needed only for 
finite hole doping x. Fictitious unit flux tubes perform 
the JW transformation. The tubes are reflected by uni- 
tary operators ul = exp[— i ^ a ■ dr\ where a is a vec- 
tor potential for tubes attached to each particleA The 
raising operator = fju^ and destroys (in the sense 
u^u] = 1) a flux tube before creating a spin particle, 
Siz = fUi - ^(1 - The physical electron = //^^ 

while = u^b^ involves a flux tube. 

The standard mean-field slave-boson method intro- 
duces unphysical states because a certain constraint is 
obeyed only on average. The present development cor- 
rects this defect. This has important physical conse- 
quences. Within the usual approach, for doping x, the ef- 
fective spin hopping matrix element ~ {J-\-2xt)/2 ~ 150 
meV which is sufficiently large to explain the strongly dis- 
persive "quasiparticle" features indicated by solid circles 
in the ARPES data and theory of Fig. Q In fact, this ma- 
trix element should be Je/2 = {J — 2xt)/2\^ where the 
minus sign reflects the well-established fact that a small 
doping X favors ferromagnetism. Near optimal doping 
Je is approximately zero. This Jg is involved principally 
in the dispersion of the so called "coherence peak," i.e., 
the "+" and "|" in Fig. ^ The "quasiparticle" feature 
reflects a fluctuation part of the hole spectrum which is 
"driven" by the constraint. The physics of the 6-particle 
hole Bose condensation is also changed when the con- 
straint is enforced. Such a condensation needs the lowest 
energy level to be renormalized to the chemical potential 
and this results in a gap equation which determines Tc 
in both the underdoped and overdoped regimes. This 
paper introduces this interpretation of the ARPES data. 
In the process some rather subtle effects associated with 
the bilayer splitting are explained. 

In analogy with the quantum Hall effect, the flux tubes 
reflected by the u] are nontrivial. That c-^ = blf^ while 



blu^ suggests that f^^ and u^^ are intimately re- 
lated. Important is the degeneracy of the vacuum. This 
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is I — however^ since H is rotationally invariant, 

the ferromagnetic state | — N/2 + 1) obtained by action 
of the total spin 5'+ on | — N/2) is equivalent. Arbitrary 
A/'/2-/-particle states are of the form iI)n/2\ ~ or 
0Ar/2-i| ~ + 1) where '0Ar/2 and (/>Ar/2-i are suitable 
field operators which create respectively N/2 and A^/2— 1, 
/ particles. Consider = h\u^^ it can be shown rigor- 
ously for large N and M ~ N/2^ f particles, 

{-N/2 + l\cj>l_,blu,i;M\-N/2) 
= {-N/2\4>l_,blfM-N/2). 

Thus the vacuum off- diagonal matrix elements involving 
are equal to vacuum diagonal matrix elements with 
this operator replaced by Here destroys a parti- 
cle, as does but the destroyed particle appears in the 
vacuum so 5'^ is unchanged. A judicious choice of inter- 
mediate states is called for. With diagonal elements is 
simply a phase factor, its particle nature only being man- 
ifest when the intermediate states introduce off-diagonal 
elements. 

Corresponding to a "flux state" , in the absence of dop- 
ing and treating the flux at the mean-field level, the 

spectrum of the / particles is, = yGe^7^^"H~"^o^^^, 
where 7^ = (cos + cos ky) while dj^ = (cos kx — cos ky). 
A T = 0, a relationship Jg = Aq = J/V2 reflects the 
unitary nature of the u^. 

Bose condensation involves the construction of a hole 
coherent state which satisfy the constraint. This is done 
by performing a SO (3) rotation. The space for a given 
site i is spanned by the basis vectors | j), | |) and 
|0). The small rotation causes |<j) W) -\- ^cri|0) while 
|0) |0) — X^cr ^cri|cr)- The rotation angle Oi constitutes 
the wave function for the condensed holes. It is necessary 
to minimize the expectation value of H with respect to 
this rotation. The result (dH/dOi) = is of the form 
of an equation of motion for hi in which the energy has 
been equated to the chemical potential /i; i.e., it amounts 
to the evident requirement that the lowest Bose level be 
renormalized to ja. The hopping term for these bosons 
is — J2Ujulbib^jUj^ and the resulting diagonal term in H 
is maximized in the negative sense when when UiOii = 
a constant, while UiO^i = ±^ where the sign alternates. 
Substituting this into — XI and making the 

same mean-field approximation as used for the J term 
results in —0'^{l/V2)^tijfjfj which is without doubt 
ferromagnetic in nature, and Je = {J — 2xt)/V2. A sim- 
ple spin pair amplitude, e.g., {flf^ J)^ is forbidden by 

the constraint, however for half filling it is well known 
that the RVB state is reflected by a flux state amplitude 
( fl ) where tt = (tt, tt). Since here the holes are sub- 
sumed into the spin particles, this result is extended to 
finite doping, i.e., with the above for 6>cri, it follows that 
the BCS pair amplitude {ctc^ ) ^ 6>2(/l/^ i.e., the 
flux state plus coherent doping equals superconductivity. 



The 0(3) rotation of H introduces new 0^ terms of 
physical importance. Since |00), reflecting near- neighbor 
hole sites, mixes with | i 0), while | ||) mixes with |0 |), 
there is a 6-particle pair term —t6'^u\ujb\b'^- + H.c. which 
results in a charge gap Kc ~ 46>^t. Rotating the physical 
electron operator produces both Oif}^ and Oiu\fl f^. 
The former coherent term implies that the condensate 
fraction is c = ^^/2 while the latter is corresponds to an 
equal number of constraint-driven incoherent fluctuating 
holes, this at zero temperature. 

In other than one dimension, the Bose levels are renor- 
malized towards ji by the —"^tijfib^b^jfj term. Using 
the usual renormalization group techniques this results 
in a gap equation. High-energy intermediate states are 
integrated out down to some cutoff K{ii). Since the in- 
termediate states involve both spin and charge degrees of 
freedom there are two contributions to K{ii) = Ks -\- Kc. 
The charge pseudogap Kc ~ 46>^t while Kg ~ ^k=7T 71/2 
reflects the width of the spin band due to both disper- 
sion and the spin BSC gap. The T = charge gap Kc 
is the lowest-energy scale and Tc ~ i^c, i-e., Tc cx c. 
Spin pairing occurs at a higher mean-field temperature 
Ts ~ J/2 a/2 associated with the opening of a pseudo- 
gap. The scenario is supported the fact that Tc (x c 
and evidence^ of a pseudogap well above Tc, in tunneling 
for overdoped Bi2Sr2Cu06+(5 and via Zeeman splitting^^i 
in Bi2Sr2CaCu08+5. With decreasing /i, when conden- 
sation first occurs, K{fi) = Kg {Kc = 0). This point is 
unstable precisely because Jg = (J — 2xt)/\/2, whence Kg 
decreases with increasing x. In the underdoped region, 
when dKs/dx < 0, the system separates into hole-rich 
regions with (J — 2xt) ~ and with hole-poor regions 
between. Doping occurs at constant ja. In the overdoped 
region dKg/dx > and implies that Kc must decrease 
and with it the condensed fraction c and Tc. Up to this 
point the holes have been accommodated within the ro- 
tated / particles. Now b particles must be added and 
this requires that /i decreases. This must increase i^r(/i) 
which slows the flow out of the condensate with increas- 
ing X. These 6-particle holes are not particularly impor- 
tant for ARPES. They are hard-core bosons which must 
avoid the existing condensed bosons. This implies a cur- 
vature to their wave functions which pushes their energy 
to at least ~ 4xt, i.e., out of range of the dispersive fea- 
tures near the Fermi surface. Their presence is seen in 
Je = {J - 2xt)/V2 since this is a mean-field result and 
all of the holes must be counted in the x. 

The physical electron Green's function describes the 
excitations of interest. Consider first a single plane. The 
operator— e.g., cj^ ^ l^il// + //^^— contains a condensate 

part proportional to \b^\ = 0/^2 and a noncondensate 
part reflected by b[. If within the Fourier transform of 
—t^cl^c- = — 2t^77:ct c- the single matrix element 
7^ is removed, then all coherent motion of holes with 
this wave-vector is suppressed and this decouples these 
two parts of Q^. If Q^^ is the Green's function calculated 
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FIG. 1: (a) Experiment"'^ (top) and theory (bottom) for the 
normal state. Theory corresponds to '^ttiQ^^ with s = 35 
meV, ^/ab — 25 meV h — 0.65, d = 54 meV and e = 35 meV 
for ky^n and k:, = Ott, 0.0257r, 0.057r, 0.0757r . . . 0.2257r in de- 
scending order, (b) The same for the superconducting state 
but with s = 6 meV, Vab — 31 meV h — 1.15, d = 50 meV 
and e = 70 meV. In both panels the strongly dispersive bond- 
ing band is indicated by a filled circle. For the normal state 
the triangles indicate the position of the weakly dispersive 
residue of the antibonding band while in the superconduc- 
tive state the | and + indicate the positions of the coherence 
peaks associated with the bonding and antibonding bands, 
respectively. The rest of the parameters are given in the text. 



without 7^, then trivially the full Q^^ = G^^ — (~2t)7^. 
A simple case would ignore both pairing and the noncon- 
densate part, whence Q^^ = {co -\-is — Jj^)/c^ leading to 

a QZ^ with poles at energy Sj^ = {J — 2ct)^^/^/2, with 
strength c, and illustrating the contribution of the con- 
densate to Jg. 

When the bilayer coupling is included, eigenstates are 
simultaneous eigenvectors of H and a reflection operator 
R; i.e., solutions, with index a = ±, are bonding or an- 
tibonding respectively. Neutron data on the insulating 
compoundsii suggest that the inter-plane J± is vertical. 
The bonding and antibonding d-wave order parameters 
must then be identical for symmetry reasons and J± is 
unimportant. With a bilayer slitting the condensate and 
noncondensate parts are decomposed by removing the 
matrix elements involving e^^ = —[2^7^ + a^{coskx — 
cos ky)'^] (rather than 7^) and the Green's function be- 
comes g^^^iu, + is) = [{G^^,^ + Gj)/{1 - e^^ [G^^,^ + Gf]). 

The full condensate contribution to Qj^^ is easily seen 



to be 



u^is-{{x-c)e^^^J-f^) 



uj-\-is-\-{xe^ 



(1) 

where included are the complications associated with both 
pairing and the bilayer splitting via e^^. As described 
above, the important noncondensate part corresponds, 
in fact, to the fluctuations of the condensate driven by 
the constraint. At this stage in the development of the 
theory this is modeled by 
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where ^mSj = (a/o;) + h{ijO — d) and where y ah ^ K. 
Then SReS/ = esgn(cc;) and d also ~ K result from shifts 
away from the Fermi surface due to the spin and charge 
gaps. In order to simulate finite temperatures, s is used 
to mimic thermal broadening, the condensed fraction c 
is treated as a parameter, and Xc is the total number 
of bosons including those in the fluctuating part of the 
condensate hut excluding the h bosons mentioned above. 
The ARPES and tunneling are simulated using the re- 
sulting g^^{uj + is). 

Normal state ARPES spectra are calculated with the 
scale ~ 25 meV with no holes condensed. With the 
other parameters x ~ 0.2, J = 45 meV, 2xct = 396 
meV and 2xct± = 79 meV, G^j^{^ + is) reproduces well 
the corresponding experimental spectrum of Feng et ali. 
Fig. nji. The honding hand lies at a binding energy of 
~ 50 — 100 meV and is quite broad. The peak which lies 
within ~ 20 meV of the Fermi surface for smaller kx is 
a remnant of the antihonding band which lies ahove the 
Fermi surface. The strength of this antibonding band 
has been artificially increased by a factor of 2 in order to 
simulate the apparent experimental fact^ that under the 
conditions with which these spectra where taken, ARPES 
couples more strongly to the antibonding states. The 
nature of the near- Fermi-surface peaks for k = [tt, 0] is 
not inconsistent with very recent data^^ for overdoped 
Bi2Sr2CaCu208. Experimentally the ratio of the normal- 
state peak and hump intensities can be varied via the 
photon energy which changes the relative coupling to the 
bonding and antibonding bands. Our simulations con- 
firm for this direction that the hump and peak are pre- 
dominantly of bonding and antibonding character, re- 
spectively. 

The superconducting state implies hole condensation 
with a condensed fraction of 0.5, and with this coher- 
ence peaks develop near the Fermi surface. The thermal 
broading reflected by 5 = 6 meV has been reduced. The 
reduction in the effect of thermal filling of the spin sec- 
tor gap is reflected by an increased K ^ 31 meV. The 
other parameters are the same (see figure caption). A 
quasiparticle honding band, again at a binding energy of 
~ 50 — 100 meV, is identified in experiment spectra and 
indicated in the theory by a solid circle. The bonding 
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FIG. 2: The decomposition of theory into bonding and anti- 
bonding bands. The sohd circle and square indicate the quasi- 
particle bands and the | and + the superconductive coherence 
peaks. The weight of the antibonding band is multiphed by 2 
to simulate experiment. The inset shows theory for tunneling. 

and antibonding coherence peaks identified by Feng et 
alX are also indicated. The decomposition into bonding 
and antibonding bands for k = [tt, 0.1257r] is shown in 
Fig. [3 The antibonding band which lies above the Fermi 
surface "pulls" its coherence peak towards the Fermi sur- 
face, leading to an apparent gap of only ~ 20 meV. On 
the other hand, the full gap ~ 36 meV is exhibited by 
the coherence peak of the bonding band. 

That the coherence peak splitting does not lead to a 
splitting of the tunneling spectrum, as it would in a sim- 
ple BSC theory, is also illustrated by the inset in Fig. [21 
This spectrum was obtained with identical parameters as 
for the body of this figure by averaging over 100 points 
in the Brillioun zone. It is interesting to note that the 
peak-dip-hump structure is reproduced with the dip at 
^ 60 meV being consistent with typical tunneling data^ 

The value of J = 45 meV reflects the gap via A = 



J/ a/2. This is small compared with J = 120 meV de- 
duced from a fit to the spin waves. However, a quantum 
correctionii reduces the latter value to J = 100 meV. 
The same correction reduces the antiferromagnetic order 
parameter by ~ 60% and assuming a similar correction 
here increases the present J to ~ 74 meV. The difference 
is accounted for by a factor of (1 — x) which arises for 
the 0(3) rotations and which was omitted in the above 
equations for clarity. 

The value of 2xct^ = 79 meV is larger that the value 
45 meV deduced directly from experiment i On the other 
hand, the larger bare value t±_ ^ 180 meV agrees well 
with the local density approximation value. That there 
is a strong renormalization of t± has been discussed 
elsewhere.^- 

While the present JW formalism is intended for all dop- 
ing levels, the present application has been to the over- 
doped regime. As described above, the present approach 
suggests that the underdoped region phase separates into 
hole-rich and hole-poor regions where the concentration 
in the hole-rich regions corresponds approximately to op- 
timal dopingii^ 
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